Spin Fluctuation Dynamics and Multiband Superconductivity in Iron Pnictides 
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Multiband superconductivity, involving resonant pair scattering between different bands, has 
emerged as a possible explanation of some of the main characteristics of the recently discovered 
iron pnictides. A key feature of such interband pairing mechanism is that it can generate or en- 
hance superconducting pairing irrespective of whether it is attractive or repulsive. The latter case 
typically leads to the superconducting gap switching its sign among different sections of the Fermi 
surface. In iron pnictides, the natural scenario is that the gap changes sign between the hole and the 
electron Fermi surfaces. However, the macroscopic symmetry of such an extended s'-wave state still 
belongs to the general s-wave category, raising the question of how to distinguish it from an ordinary 
s-wave. In such a quest, it is essential to use experimental techniques that have a momentum space 
resolution and can probe momenta of order (tTjTt), the wavevector that separates the hole and the 
electron Fermi surfaces in the Brillouin zone. Here we study experimental signatures in the spin 
fluctuation dynamics of the fully-gapped s- and s'-wave superconducting states, as well as those 
of the nodal d- and p-wave. The coupling between spin fluctuations of the incipient nearly-nested 
spin density-wave (SDW) and the Bogoliubov-deGennes quasiparticles of the superconducting state 
leads to the Landau-type damping of the former. The intrinsic structure of the superconducting 
gap leaves a distinctive signature in the form of this damping, allowing it to be used to diagnose 
the nature of iron-based superconductivity in neutron scattering and other experiments sensitive to 
spin fluctuations in momentum space. We also discuss the coexistence between superconductivity 
and SDW order. 



INTRODUCTION 

Recent discovery of a new high-temperature supercon- 
ducting family [J, [3, i, 0, i, 00, IE II has generated a 
flurry of excitement. Many of the key questions, both 
theoretical and experimental, remain unanswered. How- 
ever, it is becoming rapidly clear that the resemblance 
to the high-Tc cuprates is less straightforward and that 
a new superconducting mechanism might be at play. 

Variety of order parameters and pairing mechanisms 
has been suggested. Phonon interaction alone seems too 
weak to ex plai n high Tc. In several theoretical papers 
lol 11, H, 1^ the multiband superconductivity in dis- 
cussed as a possible explanation for the high Tc- Indeed, 
it has been long known ^14] that multiband effects can 
strongly enhance superconductivity. This is a particu- 
larly relevant in the case of iron-pnictides, which appear 
to be moderately correlated electron systems, with large 
number of Fe d-bands at and near the Fermi level. An- 
other advantage of this mechanism is that the interband 
pair interaction can enhance superconductivity irrespec- 
tive of whether it is attractive or repulsive, provided that 
the gaps in different bands have the same or opposite 
signs, respectively. The former is the familiar s-wave 
while the latter is the extended s-wave superconductiv- 
ity, or s' . This s' state, with the superconducting gap 
having the opposite sign on hole and electron sections 
of the Fermi surface (FS), emerges as a natural explana- 
tion for the superconductivity in Fe-based compounds, 
given the proximity of these materials to various nesting- 
driven spin density-wave (SDW) and related instabilities 
and absence of obvious strong attractive interaction. 



The early point-contact Andreev reflection experi- 
ments indeed indicated a fully-gapped superconductor 
[isl ITij with no indication of cuprate-like nodes, con- 
sistent with this s' picture, provided that the hole and 
electron gaps are of a similar ma gnit ude. The subsequent 
microwave |r31 and ARPES [lalSl experiments further 
fortified the case for s or s' state, also finding the fully- 
gapped superconductor with, in some instances, different 
gaps for the hole and the electron portions of the FS. In 
contrast, the NMR results are most naturally interpreted 
in terms of a d- or p-wave nodal state 20] . 

This state of affairs underscores the urgency of set- 
tling the issue of the gap structure by additional experi- 
mentation. Very recently, several papers 
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elaborated different possible experimental signatures of 
the s' superconductivity as well as of some other forms 
of superconducting order, particularly in the NMR ex- 
periments. However, in view of the complexity of these 
materials and still relatively poor quality of the samples, 
it is unlikely that a single experiment is going to settle 
this issue unequivocally. More importantly, the s' state is 
in the same symmetry class as the standard s-wave and 
there is in principle no macroscopic experiment, anal- 
ogous to the phase sensitive measurements in cuprates 
24| . which can distinguish the two in a decisive, qual- 
itative fashion. Instead, the key difference between the 
s'- and the s-wave state is in the momentum space, not 
the real space, and thus one should concentrate on exper- 
iments that have the momentum space resolution and can 
probe the wavevectors around (tt, tt), which separates the 
hole (r) and the electron (M) FS in the Brillouin zone. 

To this end, and to further expand the range of the ex- 
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perimental techniques that can be used in this regard, in 
this paper we consider the spin fluctuations dynamics in 
the superconducting state of iron pnictides. We assume 
that the system is close to a nesting-driven spin density- 
wave (SDW) instability and compare different contribu- 
tions to the damping of spin fluctuations arising from the 
Bogoliubov-deGennes (BdG) quasiparticle excitations in 
superconductors with s- and s'-wave gap. For complete- 
ness, we also consider the case with nodal p and d-wave 
symmetry of the gap function. This problem is the su- 
perconducting state analogue of the Landau damping in 
Fermi liquids. We find that the damping is qualitatively 
different in all of the above cases and thus can be used 
to diagnose the intrinsic, microscopic nature of the su- 
perconducting order. We also consider the situation in 
which the ordered SDW/AF state coexists with a super- 
conductor. This coexistence is observed in at least two 
of the compounds [1, Q and can be induced by pressure 
or doping. 



PRELIMINARIES AND THE MODEL 

The parent compound of the 1111 class of Fe-based 
superconductors has a ZrCuSiAs-type crystal structure 
[25j . with eight atoms per unit cell. The Fe atoms lie in 
a plane, same as O atoms precisely above them, in the 
adjacent rare earth (RE) oxide layer. In contrast, the RE 
and As atoms (also located above each other) are puck- 
ered out of plane in a checkerboard fashion. This puck- 
ering of As atoms is crucial for understanding the elec- 



dition, SmOi_a;Fa;FeAs and Bai_a;Ka;Fe2As2 have a siz- 
able region of coexistence of antiferromagnetism (SDW) 
and superconductivity. 



tronic structure of the compounds [ll| . It brings all Fe 
d-orbitals close to the Fermi level and creates significant 
overlap among Fe and As atomic orbitals. The result is a 
rich band structure. There are five bands crossing Fermi 
surface: two electron cylinders around M point; two hole 
cylinders plus a hole pocket around F point Ifl, 2^, 27 1. 
The 3D hole pocket is quickly filled with doping and is 
believed to be irrelevant to both antiferromagnetism and 
superconductivity. The remaining two hole and two elec- 
tron bands are almost cylindrical and exhibit significant 
degree of nesting. The natural Fe magnetism, due to the 
Hunds rule, is suppressed and one is left with a weaker 
itinerant (antiferro)magnetism, sensitive to the nesting 
features in the band structure and associated with mod- 



eratly strong correlations 

The outlines of the generic phase diagram of iron pnic- 
tides have began to emerge H, 2^ S^l : at zero and mod- 
erate doping and zero temperature there is structural dis- 
tortion and antiferromagnetic order which disappear at 
some higher temperature. This critical temperature is 
strongly suppressed by doping; at some critical doping 
structural distortion and antiferromagnetism suddenly 
give way to superconductivity. After this point, fur- 
ther increase in the doping level produces relatively small 
changes and the superconducting Tc is rather flat. In ad- 




FIG. 1: An illustration of our model: a hole band (c) centered 
at the r point and an electron band (d) centered at the M — 
(tt, tt) point interact via a short-range interaction U. 



An emerging consensus is that antiferromagnetism is 
due to the nesting between the hole and the electron 
bands. The result is the spin density-wave (SDW) formed 
by itinerant electrons. After particle-hole transforma- 
tion, the SDW instability is mathematically equivalent 
to the BCS one, caused by a logarithmic divergence 
Away from perfect nesting this divergence is replaced 
with a finite peak and sufficiently far away there is no 
instability for weak interactions. Because of proxim- 
ity to such instability, however, there are enhanced spin 
fluctuations in the system. This fluctuations couple to 
the Fermi liquid quasiparticles and can decay into an 
electron-hole pair, which leads to a Landau-type damp- 
ing. In the antiferromagnetically ordered state with lo- 
calized spins or in a perfectly nested SDW one expects 
an insulating behavior. Instead, the experiments show 
a drop in the resistivity at the SDW transition point 
followed by the metallic behavior 



31|. This obviously 



implies high degree of itinerancy and indicates that the 
Fermi surface is only partially gapped or that the Fermi 
level is located entirely outside the SDW gap. In presence 
of the SDW order, the collective excitations of the SDW 
order parameter (spin waves) can interact with the quasi- 
particles on the Fermi surface, thus inducing a Landau- 
type damping of the former. Similarly, in the vicinity of 
the SDW state, when the true long range order is absent 
but the correlation length is very large, the spin fluctu- 
ations - the incipient spin-waves of the SDW - are also 
damped by the aforementioned decay into the particle- 
hole continuum. When the system undergoes the super- 
conducting transition, this particle-hole continuum at the 
Fermi level is gapped by the superconducting order. Con- 
sequently, the new fermions - the BdG quasiparticles - 
are far less effective in damping the spin fluctuations and 
the decay rate vanishes as the temperature goes to zero. 
This decay rate of spin fluctuations carries a distinct sig- 
nature of the structure of the superconducting gap in 
momentum space and the associated BCS coherence fac- 
tors. 
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For the purposes of this paper we adopt the following 
simple yet sufRciently realistic model depicted in Fig. [1] 
We assume that the electron spectrum are described by 
two different bands: a hole (c) band at F point, and an 
electron (d) band at the M = (tt, tt) point with the same 
effective mass m. This is an approximation to be sure 
but a sensible one since the two hole and two electron 
bands at the FS of real Fe-pnictides resemble each other 
to a reasonable degree. In fact, Fe-pnictides are really 
semimetals, in the following sense: consider two bands 

ei = ed + tdCos{pxa) +tdCOs{pya) (1) 

and imagine the situation where — Ec 3> tc,td and 
we have two electrons per unit cell, mimicking the six 
d electrons of the Fe-pnictides parent compounds. The 
chemical potential is in the gap between the c (full) and 
d (empty) bands and the system is an insulator. Now, 
as the difference Sd — £c is gradually reduced while the 
electron number remains unchanged, the bottom of the 
d band at the corner M of the Brillouin zone (BZ) moves 
below the top of the c band at its center F. The electrons 
filling the top of c now migrate to the bottom of d leav- 
ing the holes in c behind, thereby creating the FS shown 
in Fig. [TJ The density of electrons filling the bottom 
of d is precisely equal to the density of holes in c band, 
hence the semimetal label for the parent compounds. Of 
course, the situation in real materials is not as simple 
as Eq. ([1]); there are four not two bands and they are 
far from being simple since their orbital content changes 
considerably as one goes around the FS pjj] . These com- 
plexities notwithstanding, the above simple picture (Fig. 
[Hand Eq. ([T])) with tc ^ td will suffice for the purposes 
of this paper. 

The relevant interband interaction is assumed to be 
the short-ranged Hubbard U . Thus, the action and the 
Hamiltonian can be written as 



First, we consider the case of perfect nesting /i = 
and compute the damping term for the spin fluctuations, 
given by the imaginary part of the electronic spin sus- 
ceptibility: 



5*0 = I drd'^r ^^c{dr - fi) c + H'' + d (dr - fi) d + H'^^ 
S,„t ^ U \ At A^r ccdd 

p- 



H'^ ^Y. '"Mp = + ^^J^ I d,d. 



2m 



uj - 



which results in: 



Imxiq + uj) oc , 



(3) 



having the familiar Landau damping form. 

Next, we consider a more realistic case and assume 
that two bands are mismatched - their centers do not 
coincide - one being shifted by /S.kp - and furthermore 
their radii differ by 2/i, as illustrated in Fig. O These 
two quantities, /S.kp and /x, parameterize deviations from 
perfect nesting within our model. The spectrum of the 
two bands is: 



2m 



Notice that the damping term is nonzero only if the two 
Fermi surfaces touch or intersect each other, otherwise it 
is impossible to excite an electron-hole pair with a low 
energy fluctuation and momentum close to M. When 
two Fermi surfaces have two different intersection points, 
the damping term is proportional to oj: 



(4) 



2m 



where Pq 



2meF - (|AfcF|/2) 




(2) 




where we have expanded c (d) band © near the top ^IG. 2: Fermi surfaces of inequivalent hole and electron 

(bottom). We can now shift the electron band to the F bands. Their centers are displaced by the vector Afcj. defined 

point, and call this "new" electron band e. One should, in the text, 
however, always keep in mind that the two bands are 

shifted relative to each other in momentum space by M = When the two bands are touching each other at a 

(tt, tt) (Fig. [T]). point on the Fermi surface, |Afci?| = ^2m(tp + /i) — 
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^^mitp — fi). For fi <^ ep but still finite, the low fre- 
quency damping term is: 



Imx{q + M, Lu) cx 



2fj, 



(5) 



SPIN FLUCTUATIONS IN A 
SUPERCONDUCTING STATE 



dHl) becomes (using Ep^^ — Ep 



epvpq cos f 



for small q) 



uj — '^p'"'^^"^^ ^ with cos 6 having to be positive for e^. 
We can now use es = 



El - Al to solve for Ejj = E: 



S{cj-{Ep+,-~Ep)) = 



1 EWE^ - A2 



vpgcost 



A2 



S E- 



A 



Once our system enters the superconducting state, the 
damping changes and becomes temperature dependent 
reflecting the opening of the superconducting gap. The 
coupling of the spin fluctuations to the new excitations 
- the BdG quasiparticles - is determined by the BCS co- 
herence factors (32). These coherence factors, however, 
have a different form for different forms of the micro- 
scopic superconducting order parameter - the BdG gap 
function - and therein lies the possibility for diagnosing 
the intrinsic nature of the superconducting state. 

Again, we first consider the case of perfect nesting be- 
tween the hole and the electron band(s). Since we are 
interested in scattering of BdG quasiparticles and in the 
processes that involve spin, the correct coherence factor 
to use is Mp^tp'+g- -I- UpTp'+g- or, explicitly: 



where 



is obviously smaller then unity; we 



Vpq cos 

will assume the condition ^ 1, as usual. The occu- 
pation number factor in Eq. /{Ep+g') — f{Ep), can 
be written as f{Ep + oj) - f{Ep) w oj^^. We then 
convert the p integration into E integration, with den- 
sity of states ^p^_i^i 1 and perform the integral over the 
(5-function. This leaves us with only the angular integra- 
tion remaining: 



(vpq cos 0) 



2E^ 



ElSP- 



(1+^ dm) 



(9) 




p p+q 

p p+q 



(6) 



where = ^(e^-)2 + {A-y and Et = ^(ef)^ + (A-)2. 

Since our focus in on spin fluctuations with momenta 
around M = (tt, tt), we are justified in dropping the in- 
traband terms. The electron susceptibility then can be 
written as 




(7) 



p p+q 



p p+q 



fm+q) - f{E^p) 



where ^^q§-^ is evaluated at E — -^==5. This integral is 
non-trivial and to perform it we have to use u/vpg <^ 1. 
Despite the ^/l — in the denominator the integrand 
is rapidly suppressed in the limit x — > 1, because of the 
exponential factor hiding in ^^q^'^ ■ It is straightforward 
to check that the first derivative at = vanishes and 
that the integrand has a maximum there. Because of 
the rapid decrease away from this point, we can estimate 
the integral with the familiar formula / ~ y(0)^/'|^r^, 

where y{cos9) is the integrand. Working to the lowest 
order in we finally obtain: 



p p+q 



\ p+q 



Et) 
P' 



^0^ 



where q <^ AI . Extracting the imaginary part gives the 
general expression for damping of spin fluctuations at 
this momentum by the BdG particle-hole excitations: 

Imx{M + q,uj) = 




X (/(iSp+g-) - fiEp))5{Lu - {Ep+^ - Ep)) . (8) 



First we consider the s' case. Because of A% ■ 



p+q 



(s' implies the repulsive pairing term) the coherence fac- 
tor is strongly suppressed EpEp+^ + ^f>^f>+q ~ ^p^p+q ~ 
^{{vpqcos9)'^{l + (-^)^), after expansion in the pow- 
ers in q (we dropped the band indices for the moment 
and put e?. = ep). The argument of the (5-function in 



/mx* (g, cx uj{vFq) 



(ef +l)2r 



(10) 



This expression gives the interband contribution to the 
damping of an incipient fluctuating spin wave in an s'- 
wave superconductor. 

It is instructive to contrast the above result with the 
standard s-wave case. After performing the analogous 
calculation for the case of a pure s-wave superconductor, 
the coherence factor for the superconductor is found to 
have a constant term ~ . Keeping only this term and 
going trough the same steps as before yields: 



ImX%q,uj) cx 



OJ 



e T 



vpq (e T + 1)2T 



(11) 



Clearly, the dynamical properties of spin fluctuations in 
s and s' superconducting phase are very different: the 
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damping and thus the decay rate of incipient spin waves 
in the s' state is substantially reduced relative to the 
standard s-wave case, by the overall factor of (note 
that q is measured relative to M) . These differences and 
the specific forms of damping should be observable in 
neutron scattering or other momentum-resolved probes 
of spin fluctuations. 

For completeness, we now consider the nodal d- and 
p-wave superconductors, another possible contenders for 
the superconducting state of iron pnictides, favored by 
the NMR experiments [IS]- In this case, the main con- 
tribution to the low temperature damping comes from 
the nodal regions. The natural excitations of the system 
in this regions are Dirac fermions with linear dispersion 
which can be seen by expanding the BdG Hamil- 
tonian near the nodes js^]. We will ignore the intrinsic 
anisotropy of the BdG-Dirac spectrum, since we do not 
expect it is to change the overall form of the damping 
term. The coherence factors for the individual nodes dif- 
fer for the intra and interband processes, but we need 
to add all the nodal contributions. The combined coher- 
ence factors are equal to lowest order (up to a numerical 
prefactor). First, we assume the PxPy c?-wave order pa- 
rameter {dxy), which means that the direction of M and 
the nodal directions are rotated by 7i7r/4 (where n is in- 
teger) with respect to each other. We expand around 
that nodes and do the momentum integral first, using 
the ^-function: 



p+q ■ 



Er. 



' de 



EdE, 



where now E ^ vp yPx + Py the linear dispersion of 
the Dirac quasiparticles. The integrand in the remaining 
angular integral is confined in the region 9 £ (7r/2, 37r/2), 
and is peaked around — n. We again use the formula 

/ ~ v(0)a/ , to estimate the lowest order contribu- 
^ 'y \y (0)1 

tion. This gives the final low temperature limit result: 



/mx'^"" (g, uj) (X uj^ql + ql 



1 



Tcosh- 



2 M 

T 



(12) 



The reason for this particular temperature dependence 
is the fact that we have restricted ourselves to the re- 
gion Lo <g; vpq- Integrating this expression over q we can 
obtain the standard result for the NMR damping rate 
l/{TiT) (X T2. 

The other possibility for a d-wave order parameter is 
a dx2-y2. In this case the nodal directions are along M 
or perpendicular to it. To the lowest order the damping 
in this case has the same form as in the dxy case, so the 
position of the nodes with respect to M is irrelevant. We 
emphasize that this is only the lowest order result, so 
in general some distinction between the dxy and c?^2„j,2 
cases is expected. 



Now let us consider a p-wave superconductor with two 
nodes along x direction. Combining the coherence fac- 
tors for the nodes gives an additional factor of sin^ 9, 
compared to the dxy case. We again use the (5-function 
to do the momentum integral first. After that the in- 
tegrand is confined in the (7r/2,37r/2) interval, but now 
goes to zero at 6' = tt. Nonetheless, we can still estimate 
the integral and the leading term behavior is the same as 
in the dxy case (with a different numerical prefactor): 



Imx^" {q, w) oc uj^ql + q\ 



1 



^Tcosh^ 



kl 



(13) 



We obtain the same result for a p-wave superconductor 
with nodes along the M direction. 

A multigap nodal superconductor can also be an ex- 
tended d and p-wave, with sign change between the gaps 
on the disconnected parts of the Fermi surface: a d' or 
p' state. The lowest order term in the coherence factor 
is then proportional to q^, as in the case of s' supercon- 
ductor, and the damping in d' and p' states is strongly 
suppressed compared to the pure d and p-waves. We see 
that momentum resolved measurements can distinguish 
the p and rf-wave superconductor from the s and s' case. 
This comes in addition to the significant differences in 
their temperature behavior. 

Perfectly nested bands generate strong spin and charge 
density- wave (SDW and CDW) instabilities. It is be- 
lieved that superconductivity appears with doping or by 
application of pressure, away from the perfect nesting, 
once the ordering in the particle-hole channel is sup- 
pressed. Thus, we now turn to the case of imperfect 
nesting. To this end, we again consider the two bands - 
a hole and an electron one, illustrated in Fig. [2l We also 
assume that the deviation from perfect nesting is small 
and Akp is of the order of q and smaller that A/vp. 
Now expanding the coherence factor in the s' case gives 
us ^{2^i + p- {q + AkF)/my{l + i^)"^). Repeating the 
same steps as before we obtain the damping: 



ImX^ {q,i^) oc uj{2fi + VF\q + Akpl) x 
[2i^+ VF\q + A^^l eT 

VF\q + AkF\ (eT-fl)2T' 



(14) 



For the s-wave superconductor coherence factor we again 
keep only term and the result is 



Imx^ (q,uj) oc 



l2ii + VF\q + AkF 



2^i + VF\q + AkF\ 

A 

e T 



VF\q + AkF\ {cT +IYT 



(15) 



Both results reduce correctly to the perfect nesting 
case {p, 0,Afci? — > 0) but generally have more com- 
plicated behavior. Increasing the mismatch of the two 
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bands drives the damping to the same form const x in 
both cases, albeit with different constants. 

In the case of d- and p-wave superconductors, the im- 
perfect nesting of this type leads to a rather comph- 
cated behavior. The point where the two Fermi surfaces 
touch (Fig. [2]) is special and its position with respect 
to the nodes determines the low-energy, low-temperature 
response. If there is a node in a close vicinity of this point 
of contact between two Fermi surfaces, the behavior of 
the system is somewhat similar to a perfectly nested sin- 
gle node case, but with M replaced by M -I- Akp. If 
this point is in one of the anti-nodal regions the response 
resembles the one found in the s-wave case. 

We pause here to emphasize again that the expressions 
for damping of spin fluctuations presented above - and 
similar other physical quantities - at a specific momen- 
tum around M are better suited to distinguish between 
different microscopic order parameters and their asso- 
ciate BdG gap functions. In contrast, the NMR relax- 
ation rate is an integral quantity and thus always includes 
both the interband and the intraband contributions. Ex- 
tracting one of them is thus a rather challenging task. 
For example, even in the s' case, the Hebel-Slichter peak 
is expected to occur [2^, due to the intraband (pure s) 
processes. 




SDW 




FIG. 3: Coexistence of spin density wave (SDW) and super- 
conductivity in our model. The chemical potential fi is above 
the SDW gap while the superconducting gap is always pinned 
to /i. 



the electron spin susceptibility % determines the decay 
rate of magnons. In the SDW-only phase we have to 
rewrite the electron operators as a superposition of the 
new eigenstates, which diagonalize the broken symmetry 
Hamiltonian. This leads to the appearance of new co- 
herence factors, reflecting the SDW order in the particle- 
hole channel [35|. Since we are interested in scattering 
of quasiparticles and in processes that involve spin, the 
correct coherence factor to use is upup^^ — v.pVj;+^: 



COEXISTENCE OF SUPERCONDUCTIVITY 
AND SDW 



(16) 



We now consider the combination of SDW and super- 
conducting orders. We will assume that the main SDW 
(particle-hole) gap opens below the Fermi level, so the 
system remains metallic, as it appears to be the case 
experimentally (Blj . This is depicted in Fig. [3] and it em- 
ulates the situation in real Fe-pnictides, where the mis- 
match among hole and electron pockets results in recon- 
struction instead of complete disappearance of the Fermi 
surface 



4- 



/)q , with 00 being 
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|34|. Then - on top of this SDW - either 
attractive or sufflciently strongly repulsive interband in- 
teraction generates at some lower temperature a new su- 
perconducting (particle-particle) gap, this time located 
precisely on the Fermi surface (see Fig. [S]). The question 
is whether spin waves dynamics can be used as a probe 
of the microscopic structure of the superconducting order 
parameter in this case as well. 

The magnon (spin- wave) dispersion relation is given by 
the inverse of the transverse magnetic susceptibility 
and is — c^q^ at the bare level, where c is the spin- 
wave velocity. Damping will be introduced by including 
magnon-electron interactions, as before. Here we neglect 
contributions from impurities and higher order magnon- 
magnon interactions; these are generally present but can 
be distinguished from the itinerant particle-hole damp- 
ing by their different and lesser sensitivity to the opening 
of the superconducting gap. Then the imaginary part of 



where ep = p^/2m—tF and ^j; = 
the SDW order parameter. 

In presence of the SDW, the unit cell doubles and M 
becomes a vector of the reciprocal lattice. We can picture 
the new band structure with two bands at the same po- 
sition in momentum space, with avoided crossing where 
the SDW gap opens. What was previously a hole-like 
(electron-like) band now has an electron-like (hole-like) 
part as well. 

One important point should be discussed here. To have 
metallic behavior and superconductivity on the top of 
SDW we have to introduce the non-zero chemical poten- 
tial for the new (gapped) band structure. This, how- 
ever, cannot be done cavalierly, since the model will not 
be self-consistent. That is the consequence of the fact 
that within the simple model adopted in this paper, a 
gap on the Fermi surface is typically energetically prefer- 
able to a gap below or above. To remedy this and sta- 
bilize the SDW gap below the Fermi surface, we have 
to include new terms in our Hamiltonian, reflecting de- 
viations from idealized bands and the lattice effects in 
real iron pnictides. This, unfortunately, comes at the 
cost of significantly more complicated calculations and 
an entirely obscured physical picture. For the purposes 
of this paper, this is too great a cost and we avoid it 
by using the following approach: we calculate the SDW 



wave coherence factors at zero chemical potential, when 
Fermi surface is completely gapped. Next, in order to 
account for the metallic behavior, we reintroduce a small 
Fermi surface, which now can undergo a superconduct- 
ing instability. The physics behind this approach is clear: 
the SDW coherence factors determine the vertices that 
couple the magnons and the BdG quasiparticles. The 
magnons are the Goldstone modes of the SDW and this 
coupling should be of a gradient type, because of the bro- 
ken symmetry of the ordered state - the long wavelength 
twist in the spin direction of the SDW should come at no 
cost in energy. This remains true independently of the 
position of the SDW gap. Adopting this simplification, 
the SDW coherence factor is just ~ {vFqcosOy . 

Starting with the case without superconductivity, the 
damping is: 



For the case of a pure s-wave superconductor we get 



Imxiq, ^) °^ 



Vpq 4(/i2 _ 02)2 



(17) 



In the case of coexistence we have to calculate anew the 
superconducting coherence factor 



1 



^p-^p+q 



(18) 



In the vicinity of the SDW gap, the new quasiparticles are 
almost an equal mixture of the two initial bands, c and d. 
Turning on the chemical potential moves the Fermi level 
to the upper band; as discussed above this is to be done 
advisedly. If the chemical potential is not very small, 
the content of the quasiparticles at the new Fermi level 
will be almost exclusively from only one of the old bands 
(either purely an electron-like or a hole- like). With this 
assumption we calculate the superconducting coherence 
factors. 

Again, we first consider the s' case. Working to lowest 
order in lj / q we get: 



Imx{q,oj) « uj{vFqf 



e ^ 



(ef +l)2r 



(19) 



This expression gives the interband contribution to the 
spin wave damping in an s'-wave superconductor. How- 
ever, here we cannot neglect the intraband contribution 
as we did earlier, since now the momenta M and are 
effectively equivalent by virtue of the umklapp scattering 
off the underlying SDW modulation - consequently both 
inter and intraband terms are contributing. Calculating 
the intraband term gives: 



Imx{q,uj) « uj{vFq) 



A 

e T 



(eT + 1)2T 



(20) 



and we see that the this is in fact the leading term in the 
limit g — !■ 0, masking the contribution of the interband 
damping. 



Imx{q,uj) « uj{vFq)- 



A 

e T 



(ef + 1)2T 
The intraband term is of the same order. 



(21) 
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s' SC 
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perfect 
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Loqe T 
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nesting 
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T cosh2 M 


T cosli2 M 














imperfect 


LJ 








nesting 




V 









TABLE I: The summary of our results for damping of spin 
fluctuations by particle-hole excitations in several different 
superconducting states. Here a, = VF\q + Akpl ■ 

Thus, in the region of coexistence, the dynamical prop- 
erties of SDW magnons in an s and s' superconductor 
are difficult to distinguish. The damping at momentum 
around M is due to both inter and intraband scattering 
and since the latter are, of course, unaffected by the rel- 
ative sign of the different gaps, it would be difficult to 
observe any significant difference in experiments. 

Table U summarizes some of our results for the reader's 
convenience. 



A SIMPLE MODEL OF MULTIBAND 
SUPERCONDUCTIVITY 

A significant part of our focus in this paper was on 
an s'-wave superconductor, in which the Fermi surface of 
Fig. [1] is fully gapped but the gap function in its hole 
pocket has the opposite sign relative to the one in the 
electron pocket. Furthermore, our paper is clearly de- 
signed to arouse interest in the experimental community. 
Consequently, for the reader's benefit, we consider here 
a rather basic picture of s' multiband superconductivity 
and some of its features, appropriate for our simplified 
model of Fe-pnictides (for a more theoretically inclined 
discussion, 



the reader is referred to [34|)- With appropri- 
ate modifications, the same considerations can be easily 
adapted to the d' and p' cases. 

First, the model ^ discussed so far is not sufficient. 
The interband interaction Uccdd will induce the SDW 
and promote strong spin fluctuations but will not by it- 
self lead to superconductivity from purely electronic in- 
teractions [131 . For that possibility to enter into play, one 
must consider the interband pair resonance interactions 
of the type Jccdd+h.c, a Josephson-like term in the mo- 
mentum space which scatters pairs of electrons between 

mail. 



and d bands 



Such terms are typically present 
in multiband systems and their size in Fe-pnictides is sig- 
nificant 34 1 . 
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The key feature of this interband interaction J is that 
it can drive the system superconducting - or enhance 
the already present intraband superconductivity - irre- 
spective of whether it is attractive or repulsive. In the 
former case, the intraband gaps will have the same rela- 
tive sign while in the latter, which is probably where the 
Fe-pnictides belong, this sign will be different. This is 
the origin of the s' superconducting state. 

In the weak-coupling theory, however, the interband 
repulsive interaction J can drive the system supercon- 
ducting only if it is large enough, i.e., > U1U2, where 
Ui and U2 are the repulsive intraband interactions in the 
hole (c) and electron (d) bands. Such a sizeable J is not 
something that is easily found, since, generically, J is sig- 
nificantly smaller than the intraband Coulomb repulsion 
Ui, U2; indeed, this also appears to be the case in Fe- 
pnictides [34]. Here we suggest a mechanism that could 
potentially solve this problem. The interactions that en- 
ter the condition for superconductivity are not the bare 
Coulomb terms but some appropriately screened interac- 
tions. The screening length is kd = 47re^no(g), where 
no(<?) is a polarization bubble. The momentum region 
for the J pairing term is around AI and the main con- 
tribution to the polarization is a mixed bubble. On the 
other hand the main region for the Ui^2 terms is around 
and the screening is due mostly to the usual single-band 
bubbles. K£)(M) has a rather dramatic evolution with 
doping [ll[ (in contrast with K]j{0)), and goes to for 
high doping levels. Let's now confine ourselves to the 
region in which the interactions are frequency indepen- 
dent. This limit is constrained by the time r it takes an 
electron to traverse the Debye screening length r « ^ 
where — = 4^eWo ' interaction is frequency 
independent for t^; <C ^. Outside this region it has com- 
plicated frequency dependence and we can set J = (as it 
is done in the extensions of the BCS theory which include 
Coulomb interaction, where J = const for lj < loc and 
J = for Lo > ivc, where wc is the Coulomb frequency 
cut-off). For k_d(0) > kd{M), the frequency range of J 
is smaller and this could lead to logarithmic suppression 
of i7i_2 akin to the reduction in the Coulomb pseudopo- 
tential caused by phonon retardation effects within the 
standard BCS-Eliashberg theory [sl] . 

We now illustrate the above arguments with an ex- 
plicit calculation. In multiband superconductivity there 
are at least two gaps which vanish at the same critical 
temperature (if there is a non-vanishing pair resonance 
term J). We can write down the equations determining 
the zero temperature gaps and the critical temperature 
[fl [37! . To study the effects of the retardation we divide 
the energy range into two intervals: e G (0,a;ci) and 
e € ((jjci, ti^C2), where ujci and 0Jc2 are the frequency 
cut-offs for J and Ui^2 respectively. For each band we 
define two gaps, A^i and Ai2 (i = 1, 2 is the band index), 
corresponding to these two regions. Assuming parabolic 



bands as before, the gap equations are: 

tanh(^£;ii/2) 



All = - Aii^ 



de- 



Jo ^21 

tanh(^^i2/2) 
Ai2Afi / de 



E12 



tanh(/3gii/2) ^ 

A12 = - Aii^i / de- 

Jo -^11 



A„„/ "-"f"/^' d. (22) 



where Eki — y/f-li + A^^. To simplify the notation we 
have introduced new variables /Xi_2 = A^i, 2^^1,2 and Ai^2 = 
Ni^2J, where Ni is the density-of-states (DOS) in band 
i. Both /ii and Ai are positive corresponding to repulsive 
interactions. The above equations are coupled with the 
analogous ones for A21 and A22. The interactions J, Ui, 
and U2 appearing here are the angular averages over the 
Fermi surface [34| . 

To proceed, we concentrate on the condition for a non- 
trivial solution in terms of A^^. This gives an equation 
for the integrals Ii and I2 that appear in Eq. ((22|) : 

^ ^ tanh(/3£;ii/2) ^^^ 
Jo 



En 

^ tanh(/?gi2/2) ^^ 



(23) 



These integrals are positive-definite and close to the 
critical temperature can be approximated by Ii — 
log (1.13ci;c'i/rc) and I2 = log (i^C2/<^ci)- If we now mo- 
mentarily suspended the above "retardation" effect of J 
relative to Ui^2 = ^ci) and take the cut-off of the 

energy integration to be the same everywhere /2 ^ and 
h becomes 



r (mi + M2) + y/il^i - /i2)2 + 4A1A2 , . 

2(AiA2 - ^Ilfi2) 

The physical solution exists only if the right-hand side 
is positive, which translates into A1A2 > /^iA*2. Solving 
for the critical temperature we get: 



_ (mi +1^2 ) + \/ (l^l -M2)^+^^l-^ 
Tc ^ {l.l3uJc)e 2(AiA2-MlM2) . (25) 



The full expression for Tc, obtained from Eq. 
with UJC2 > ^ci and under the same conditions /i > 
and /2 > 0, is more complicated 



T, = (1.13wci) X 



(26) 



(Ml+M2+2J"2MlM2) + y(Ml-M2)^+4AiA2(l + /2Ml)^(l + /2M"2)^ 
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As it turns out, however, once the retardation effect is 
restored, all the reader needs to do to obtain a good 
approximation to ([^5]) is to simply replace in (|25p /ii_2 — > 
fii 2, where 2 ^-^^ given by: 



Ail, 2 



l + Ail,2log(cJC2/t^Cl) 



(27) 



Therefore, now A1A2 only has to be larger than 
or, in the original notation J > yyUfU^, to make s' 
superconductivity possible. Notice that, in a similar vain, 
any reduction in Ui,2 arising from the interband phonon 
attraction would help superconductivity as well. 

One of the consequences of Eqs. (|22p is that the overall 
magnitudes of the hole and electron gaps, |Ai| and IA2I, 
are generically different, as soon as the c and d band 
parameters are not the same. Several point-contact An- 
dreev reflection (PCAR) measurements 1^ lH, however, 
show a single gap. Obviously, one explanation is that 
there are two gaps, but they have similar magnitudes in 
this particular system and cannot be resolved in a PCAR 
experiment. 




0.02 0.04 0.06 0.08 

(N2-N1) /N2 

FIG. 4: The illustration of the relatively mild sensitivity of 
the two superconducting gaps to differences in the DOSs be- 
tween electron and hole bands {{N2~ N2)/N2). The gaps were 
computed from Eq. (|22|l . without retardation effects, using 
fii — 0.3 and Ai — 0.5, for up to 10% difference in DOSs, and 
are measured in units of W~^ujc- 

In contrast, the ARPES experiments do in fact show 
two different gaps in some of the other Fe-pnictide super- 
conductors ISj, [l9| . Note that the main role of different 
band parameters arises through their different densities- 
of-states (DOSs), the parameter which is independent 
of doping for 2D parabolic bands. Naturally, the real 
hole and electron bands in Fe-pnictides are neither ideal 
parabolas nor identical to each other (nor entirely 2D for 
that matter). To study the effects of different DOSs for 
the two bands we solve numerically Eq. (|22p (without 
the retardation effects) at zero temperature for differ- 
ences in DOS of the two bands of up to 10%. The result 
is shown in Fig 2] and we see that the hole and electron 
gaps diverge away from one another rather slowly. This 
illustrates the relative lack of sensitivity of the gap func- 
tions to differences in the band parameters, retroactively 



provides another justification for our assuming that they 
are the same, and, within the framework of the s'-wave 
state, can explain the experiments that seemingly point 
to a single-gap superconductivity. 



DISCUSSION AND CONCLUSIONS 

In this paper we investigated the dynamical proper- 
ties of spin fluctuations in superconductors with vari- 
ous forms of the microscopic order parameter, both in 
the vicinity of and within the itinerant SDW state. Our 
main emphasis was on finding ways to distinguish differ- 
ent types of superconducting order by concentrating on 
the features in momentum space, since one of the leading 
candidates, the extended s'-wave superconductor, does 
not differ from an ordinary s-wave by the overall symme- 
try. We demonstrated by explicit calculations that the 
momentum as well as the frequency dependence of the 
spin fluctuation decay rate can be used to distinguish 
among different states of a multiband, multigap super- 
conducting system, even in the cases when the tempera- 
ture dependence is similar; some of our results are sum- 
marized in Table m Among other applications, we expect 
our findings to be particularly useful in the ongoing ef- 
forts to establish the symmetry of the order parameter 
in the iron-based high-temperature superconductors. 
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